General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA TECHNICAL MEMORANDUM 


NASA TM-77804 


ALGORITIIM FOR SPACE-TIME ANALYSIS OF DATA ON GEOMAGNETIC FIELD 


N. V. Kulanin 


(RAiJA-TM-77U04) ALGOBlIHH FOR SPACE-TIiiE N85-176C9 

ANALYSIS OF DATA CN GEOfiAGNEIJC FIELD 
(National Aeronautics and Space 

Administration) 27 p HC A03/EF A01 CSCL 12A Unclas 

G3/64 13719 


Translation of; "Algoritlim pros transtvenno-Vremenno Anal iza 
Dannykl'i o Geomagnitnam Pole,” IN; "Algor itlimy i result 'taty 
obrabotki dannylch v MTsD,” (Algor j tints and Results of Data 
Interpretation in IffiC) , edited by V.P. Golovkov and Yu. S. 
Tyupkin, Soviet Geo-physical committee, USSR Academy of 
Sciences, Moscow, 1978, pp 1-5, 44-54 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
(WASHINGTON, D.C. 20546 NOVEMBER 1984 


JTAHDARD TIUE PACE 


1. R^pott Nn. 

NASA TM-77804 

Government Acceiilon No» 

J. R^cJplenl** Cololop No, 

Ti»l* and Subllll* 

ALGORITOM FOR SPACE-TIME ANALYSIS OF DATA ON 
GEO^iAG^^ETIC FIELD 


6, Performing Orgonliotlon Coif* 

7. AutSor(«) 

N.V. Kulanin 

8. Parformlng Of goni i ollon Report No, 

10. WoA Unit Ho. 

1). Contfoci Of Gront No, 

NASlV-4005 

9, P cfft-triinp Oiyoni lotion Add^n* 

Leo kanner Associates 
Redwood City, California 94063 

13, T,p. cf Rrpori o^<^ P*i!p^ CoArf.J 

Translation 

17. Sponioiing Agency' Momt and Addies» 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

1^. Sponrofing Agency Cod# 

15. Skipptrmcniary Notes 

Translation ofj "AlgoritJim prostranstvenno-Vremenno Analiza Dannykh o 
^oma^itn^ Pole, "IN: "Algorithy i rezult'taty obrabotki dannykh v 
CAj^oritlims and Results of Data* Int^^^^ in ilVDC), edited 

y Y*P* Golovkov andYu. S. Tyupkin, Soviet Geo-physical ce mmittee 
Academy of Sciences, Moscow, 1978, pp 1-5, 44-54 USSR 

16. Abstioct 

The algorithm described, calculated for the execution of the space-time 
analysis of large volumes of diverse data quite arbitrarily distributed 
on the surface of tlie Earth and in time and realized in the form of a .. 
program for the BESM- 6 computer, is now at the stage of completion of the 
verification of its applicability for analyzing a large part of the data 
contained in the Veynberg catalog with the involvement of some independent 
resuD.ts pertaining both to epoclis earlier than those contained in this 
catalog and later epochs, in the 20th century. 

1 

17. Key Wo)d» (Selected by Author(O) 

IS. Dl ilti bulior^ Stolcmont 

Unclassified - Unlimited 

19. SrCUfIty Closstf. (of tbis fcport) 

Unclassified 

70. SnCijrjfy Cictftif. (of ifni pogo) 

Unclassif i ed 

71* No. of Poge» 72, PrSc« 

18 




•.lA Cfi.UC 








ALGORITHMS AND RESULTS OP DATA INTERPRETATION IN WDC 


V.P. Golovkov and Yu.S. Tyupkln 
Soviet Geophysical Committee 

FOREWORD 

The broad use of computers in geophysics opened up two possi- 
bilities for researchers: construct complex theoretical models, 
on the one hand, and interpret large volumes of data, on the other. 
The first of these possibilities outpaced the other in realization 
because data reduction and data input into computers has thus far 
remained the bottleneck in the overall system, where the instru- 
mentation has the design-wise obsolescent analog data output. 

This collection of articles contains the results of studies in 
both the methodological and interpretational areas in one way or 
another associated with the interpretation of large volumes of 
geophysical data of the most diverse specialization. Granted all 
the narrowness of the problems solved in each study, as a whole 
they illustrate to some extent the possibilities that appeared 
before researchers when both large volumes of data in machine- 
readable form and computers with sufficiently developed standard 
or specialized mathematical software were made available to them. 

Naturally, these possibilities have surfaced first of all in 
the World Data Center; firstly, it has computers, secondly, a data- 
base, and thirdly, since it is the data processing laboratory of 
the CAPG [expansion unknown], it engages in working out methods, 
algorithms, and packages of specialized geophysical programs. 

In preparing this collection, we hoped that it will better 
display the possibilities of the World Data Center than the center 
^Numbers in the margin indicate pagination in the foreign text . 



catalogs and advertising brochures. Additionally, all the data 
used in the studies, just like the packages of programs, are 
accessible to broad classes of users at the World Data Center in 
the USSR and abroad. 


V.P. Golovkov 

ALGORITHM FOR SPACE-TIME ANALYSIS OP DATA ON GEOMAGNETIC 

FIELD 

N.V. Kulanln 

Several algorithms for space-time analyses of the geomagnetic 
field have been published in the past decade, including some using 
a system of basis functions differing from a system of spherical 
harmonics [1], as well as algorithms of spherical harmonic analyses 
admitting a sequence of spherical harmonics distinct from the canon 
leal sequence [2,3]. 

Given the unquestioned advantages of each of the above-men- 
tioned algorithms, for the most part they impose on the original 
data requirements that become burdensome when a conversion is made 
to temporal integrals encompassing the ante-Gausslan epochs, where, 
as a rule, data of a different type (.2), 1,7* ) were quite arbitrarily 
(and often very unsuitably) distributed in time and space. 

This study proposes an algorithm for space-time analysis im- 
posing virtually no special requirements either on the composition 
of the initial data, which may be quite diverse ,T, 

S. , ^ % ) or on their distribution in time and space. The 

primary constraints figuring in the specific realization of the 
algorithm on a computer stem exclusively from the limited possi- 
bilities of the computer involved. 

The algorithm uses and to a certain extent develops the ideas 
advanced in the studies [4,5]* The main "traditional" aspect of 



the algorithm can be regarded as the use in the system of basis 
functions of a canonical sequence of spherical harmonics (here it 
is provided, as a variant, that some groups of coefficients and 
their time derivatives be recorded with a certain orientation 
toward the Jukutake method [6]). 

Description of method. Underlying the algorithm is the re- 
presentation of the geomagnetic field potential on the spherical 
surface of the Earth with radius R; at the point {d^A) at the 
instant of time ^ (by limitation to Internal sources) in the 
form of a sum of the series [7] 



supplemented by the representation of the secular variation of 
the coefficients in the neighborhood of the central epoch 



wher e k (i) 

is the S-th derivative in time at the point t = tQ. 

Based on Eqs . (1) and (2), the expression for the components 
of the geomagnetic field intensity can be written in the form 

( 3 ) 

where f—X— , .V— or ^ is the intensity component, 

is the corresponding basis space-time function, 
the k-th in the canonical sequence. 




here the subscript is associated with/^* ^ and ^ as follows: 


for 


/X - I , 

m B 0, 

a 0 

H = I 


I 

0 

I 

2 


I 

0 


I 


I 

I 

0 

SfTL + /J i 





O 

I 

) 

I 

I 

I 

^ Hm, + 3 


i ' 

I ' 


2^+ 2 


I 

I 

0 

2 + 3 p/ 



^ * » « » « « 


(for ) 

) 

2 

0 

' • • » 1 t 

0 

f -f' ' 



and » tifr^ ( + 2 ) ( S,n + I ) , 

Thus, we turn to the linear analysis, in which it is re- 
quired to determine the unknown parameters representing 

the coefficients §n* * and their time derivatives 

(both at the instant of time ) . 


In order to apply Eq. (3), all the initial data must be re- 
duced to the form of the components X , J:/ or Z in accord- 
ance with the type of data; this is done by iteration, in which 
the (n + l)-th step is determined by the formulas 

JlL V 

T_ 

r/M^ r' 

i\' u' t' T' 

where Jj »n > J. * ' are the results of measuring Z and 
T , respectively, and JJ are the components 

computed at the n-th iteration step, ///'">- > 

-V (XS‘'^y-+ ; here *' are specified on the 

basis of (3) in accordance with the given initial values of the 
coefficients and their time derivatives at the Instant 

t~to 



We note that the use by T' of the fourth or fifth of the 
transforms (^1) is determined by the value of the polar angled. 
In place of (3)j we now can write 


where is the result of iteration (left-hand side of (^)) at 
the ^^th of the total number of measurement points , where 


. We note that the functions (which we 

assume, just like the variable Ee. , to have already been multiplied 
by the square roots of the measurement weightings formed on the 
basis of Eqs. (4)) generally speaking are nonorthogonal in the 
given domain (and all the more so on the set of points (Tl)* 


Eq. (5) is approximate even for a fairly large , since the 
quantities are subjected to the influence of measurement errors. 
Based on these data we can only approximately determine the coeffi- 
cients c{h . 

One of the methods of determination is to solve the system of 
conditional equations 




( 6 ) 


under the condition 

(least-squares method) . The coefficients thus obtained con- 
tain the errors 4 (where are the true values of the 

coefficients that we introduced from considerations of convenience 
in manipulations (see also below Eq. (15)), whose value generally 
depends on the choice of the number of terms in the series in 
Eq. (6) — of the quantity fi . 


The question suggests itself: how can one choose the optimal 
value of /I and further, how can one refine the coefficients 



within the limits of the selected spectrum such that = 7 

This problem is discussed below with reference to the concepts In 
the study [5]. 


Let us turn from the system of functions to the func- 

tions (yj that are orthogonal on the given set of points 

!%] . i.e., that satisfy the condition 




(r) 


Orthogonallzatlon is carried out according to the Gram-Schmidt 
formulas : 




( 8 ) 




The coefficients of orthogonallzatlon providing for Eq. (7) are 
equal to 



(^P) ^p) 


(9) 


Single application of Eqs . (8) and (9) when solving large- 
scale problems ( ft and V are large) may not provide for the 
rigorous orthogonality of the function owing to the accumu- 

lation of truncation errors. 


Here we propose an iterative process of orthogonallzatlon 
differing somewhat from the process given in [5], since the 
latter, as can be shown, requires additional correction. 

The next step in the iterative process is determined by the 
formulas : 
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( 10 ) 


(Pj 

c'H. P V^<^vV 

We can show that the iterative process can lead to the at- 
tainment of rigorous orthogonality, namely as 


/HQ 


cf) ob j ^ d ^ y 

But the question of the size of the domain of its convergence can, 
most likely, be solved only empirically on the basis of trial 
calculations for each specific case. We note only that, in addi- 
tion to the excessive value of M and a/ , the size of the domain 
of convergence is unfavorably affected by the poor causality of 
the system of normal equations, which can be obtained on the basis 
of system (6) . 


We can now assume that we have a system of functions (^) 
that are orthogonal on the given set of points l^e] . Here 
Eqs. (6) take on the form 


z c/,% 


( 11 ) 


By expressing in (11) the functions *^a in terms of G-n in 
accordance with (8) and using the linear independence of , we 
get 



(12) 
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where 




( 13 ) 


We note that the values of J^ux do not depend on the choice 
of /t , By multipylns (11) scalarly by , we find ( £' is 
the vector / ) : 


a 




ilM) 


The values of CZ^ are also not dependent on the choice of /^9 

^ , although system (11) is Inconsistent and the equality between 
£c and the sums is approximate.^ 

n \ ( (T 

In accordance with (3)j let us determine the exact (true) 
values of the coefficients as satisfying the system of equa~ 
tions when there is a hypothetical expansion of the Initial system 
of points fi + fv 'io . as a/'— ~ 

,fo4 Cfc)‘^e-^e (15) 

where Ac. is the error of the quantity £& due to an error of 
measurement at the point '^c . 


Based on data of N measurements, the following system of equa- 
tions can be set up 

I ( 16 ) 

Y ~ t 

’ The more or less weak dependence of JSu( and iX-/r on fc can occur 
as a consequence of the nonlinearity of the expression for the 
right-hand sides of the conditional equations (6) in terms of the 
functions (see (^)). 



By subtracting (16) from (15) » we get 

where a c£,i'^d!( -c/n , ^ 

On the other hand, from (16) It follows that 


therefore 


£4 < <#* (icYL-i 


(17) 


( 18 ) 


and 


K=y,.„.,V (X9) 

Prom ( 16 ) and (17) and referring to (8) and (13) we find 

where the quantities are determined by Eq. (19)* 

In most applications and, in any case, when the space-time 
algorithm is applied to data from the ance-Gaussian epochs 




( 21 ) 


for the vast majority of measurement stations. 


/50 


Referring to (21) in Eqs. (I 8 ) and (19) and admitting only 
a weak correlation between etc G-i C^e) and , we find 


c5a' " ^ 




( 22 ) 
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In any case, for , 

where d and 6a'„ are the root-mean- square errors of the quantities 
and , We note that S does not depend on ^ , since the 

weightings of the measurements were allowed for, as assumed, in the 
quantities [symbols not given in text] . 


Let us elucidate how the errors of the coefficients at the 
functions change if we convert from system (16) to the 
system of conditional equations 



(23) 


According to (12) and ^20) 

Hence, from (20) it then follows that 


or 




(24) 


where , Prom (7):> (19) ^ and (21) it follows that 




(25) 


where is the mathematical expectation of ,x. , and (24) in 

place of (25) leads to 




z oz 

Q.f 



( 26 ) 


Prom this it follows that in order for the quantity 
to become smaller in the transition from fl to H + x , it is nec- 
essary and sufficient that the inequality 
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( 27 ) 




be satisfied. 


In the first approximation we assume that the quantity /51 

f/ 

when there is a change in /t undergoes random 
fluctuations about the mean value equal to zero. Then the re- 
quirement of a decrease, on the average, in the quantity 

Mifh d* )'*7 when there is an increase in leads to the 
condition 


CL K.*t 



( 28 ) 


We will assume the quantity (t/t#/ to be ''random” with respect 
to the fixed quantity > l*e., we assume that 

Alfa,,, -J , CD fa**y7- 

Thus, Eqs . (19) and (21) allow us to consider the "distribu- 
tion” of the quantity to be norm (here we assume that // is 

large enough) i from (28) it follows that 

where {> is the minimum admissible probability of inequality 
(28) being satisfied, and ^ is determined from the equation 

where is an integral error function. 
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Thus, for example, S (0,93) - 0,4, S (0,35) . o,5. 

In contrast to the study [5], by solving the problem of mini- 
ml25ation , we provide for the second approximation. 


Here we assume 




where fl„ is the largest value of for which inequality (28) 

is still valid. We also introduce the notation: 




CCm: 








o-'au f 


f r.ir\ f 


CL 


<*■/ 


Then inequality (27) can be written in the form 




(29) 


By solving the corresponding quadratic equation, we determine 

(30) 

Analogously to the solution of inequality (28), we assume /52 

the quantities Z' and to be fixed, and [symbol not given] 

and [symbol not given] to be random, such that 
* 


Thus, the solution of inequality (29) can be reduced to esti- 
mating the probability of the random variable Z entering the do- 
main outside the random interval ( -</ ). We note that Z and 

— as linear functions of the coefficients <ZJ ^ — are normally 
distributed and at the same time are uncorrelated, since they 


12 



contain the coefficients with differenb subscripts, leading to 
the conclusion that 1 and $ are Independent. By applying the 
theorems of addition and multiplication of probabilities [8] and 
with reference to (30) and (3l)j let us find the desired probabi- 
lity in the form 


5_c'//7 ^ 


As a result., by referring to (12) and (25) j we get the values 
of the coefficients refined in the sense of an approximation 
to : 






(33) 


[subscript in left-hand side reads "refirjed"] 


and also 


where 




when /K ^ or when 
in the remaining cases 


Here is the boundary value of the subscript ff determining 
the transition to the optimization regime defined by Eqs . (32), (33) j 
is the minimum probability that inequality ( 29 ) will be satis- 
fied. Since the refinement of ceV was carried out successively 

times (but not more than times), l.e., repeatedly, 

it is best to set ^ 0,5. 

We note that by forming with the aid of the algorithm, 

by using Eqs. (12) and (13) we actually replace with the 

Jit in Eq. (12), which means that in the transition across the 
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value of i (Eq, (33) from ^6< « d a replacement in Eq. (12) /53 

of the form f^cC^ , which ultimately requires 

the substitution in Eq. (32) of the form S*-’' €Li 

The possibilities of the procedure of refining the coeffi- 
cients cIk determined by Eqs. (32) and (33) are now being inves- 
tigated with test models. 

Thus, on a model utilizing the quasiuniform filling of the 
surface of a sphere with "measurement stations" with the super- 
imposition on the data synthesized at these stations, pertaining 
to the same epoch (1965), of the error typical of the measurements 
in the l6th through l8th centuries, we obtained - 25 when 
= 35. 

Let us briefly describe the order of execution of the algo- 
rithm: orthogonal! zat ion (Including repeated) j iteration in accor- 
dance v/ith Eqs. (4)i computation of the coefficients, cl^ and of 
their errors; determination of fin*. ; when 7R. (where a is the 
last overestimation of the possible length of the series (3))j 
expansion of the list of basis functions and repetition of the 
preceding steps; selection of the initial data in accordance with 
the ratio of nonclosure to the root-mean- square error; and repeti- 
tion of the preceding steps (singly); and finding u K yfOV#, J • 

The algorithm described, calculated for the execution of the 
space-time analysis of large volumes of diverse data quite arbi- 
trarily distributed on the surface of the Earth and in time and 
realized in the form of a program for the BESM-6 computer, is now 
at the stage of completion of the verification of its applicabi- 
lity for analyzing a large part of the data contained in the Veyn- 
berg catalog [9] with the involvement of some independent results 
pertaining both to epochs earlle~- than those contained in this 
catalog and later epochs, in the u ' century. 
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In addition, the algorithm has been .applied to an analysis of 
paleomagnetic data of the Bryunes epoch [10]. 

In conclusion, the author expresses his gratitude to N.P. 
Ben'kova, V.P. Golovkov, and T.N, Cherevko for valuable comments 
and discussion. 

The author is also indebtedness to S.I. Braginskiy, under 
whose supervision he began developing the algorithm described 
above . 


15 


REFERENCES 


1. Pushkov, A.N., E.B. Paynberg, M.V, Plskin, and T.A. Chernova, 

"Space-time analysis of main geomagnetic field by method of 
decomposition into natural components," In the collection: 
Prostranstvenno-vremennaya struktura geomagnitnogo polya 
LSpace-Time Structure of Geomagnetic PleldJ, Moscow, Nauka, 
1976. 

2. Bazar zhanov, A.D. and G.I. Kolomiytseva, "Improvement in the 

analytic representation of the secular variation," Geomag- 
netlzm 1 aeronomiya 7/5 > 868 ( 1967 ). 

3 . Kolomiytseva, G.I. and A.N. Pushkov, "Analytic model of main 

geomagnetic field for an Interval of 2000 years," in the 
collection: Prostranstvenno-vremennaya struktura geomagnit- 
nogo polya , Moscow, Nauka, 1976. 

4. Cain, J.C., W.E. Daniels, S.J. Hendricks, and D.C. Jensen, 

"An evaluation of the main geomagnetic field, 19^0-1962," 

J. Geophys . Res. 70/15» ( 1965 ). 

5 . Mishin, V.M. and A.D. Bazarzhanov, "Selection of spectrum of 

Legendre polynomials approximating the observed [text miss- 
ing] -field," in the collection: Geomagnitnyye Issledovanlya 
[Geomagnetic Studies], Moscow, Nauka, No. 8, 1966. 

6. Jukutake, T., "Spherical harmonic analysis of the Earth's mag- 

netic field for the 17-th and l8-th centuries," J. of Geom. 
and Geoelectr. 23/ 1, 11 (1971). 

7 . Yanovskiy, B.M., Zemnoy magnetism [Terrestrial Magnetism], 

Leningrad, Izdatel'stvo Leningrad State University, 

1 , 86 ( 196 ^). 

8. Ventsel', Ye.S., Teorlya veroyatnostey [Probability Theory], 

Moscow, Nauka, p. 37 j 1969. 

9 . Veynberg, B.P. and V.P. Shibayev, "Catalog of results of mag- 

netic determinations on the globe from 1520 to 19^0," Ser- 
lya rukopisnykh trudov In-ta zemnogo magnetizma [Series 
of Manuscript Works of the Institute of Terrestrial Magne- 
tism], Leningrad, 19^14 Issue 5> Vols. 1-3 (with rights of 
manuscript) . 

10. Ben'kova, N.P., V.P. Golovkov, N.V. Kulanin, and T.N. Cherevko, 

"Analytic models of the geomagnetic field in the Bryunes 
epoch," in present collection. 


16 




ci 




hi 


am 


er^S 

p. 5 

ji 


I 


‘ 

o 

Q 

» 

f 
• •» 
•O 

’ 4 > 

V 

'"S 




M 

o 




M 

O 

& 

V 

■§ 

5 


fi 


U 

o 


f- 

U * 

m 

4» 

() 


I 

« 

!■« 


•■3 

!a 

:S 

<»■ 


r-t^ 




« I 

5 


|5 

! « 

r k3 

'9 

I « 

(I 

V*H 

r 

«• o 
f* w 


d 

o 

3 

4 * 

41 

u 

S' 


s 

Q 

M 

O 


I 

f I 



I 

I 


i 


V V 


ORIGINAL* PAGE IS 
OE POOR QUALITY 



{ 


CD 

o >: 

:si rq 

p n 

o t? 

3 ■& 


» a: o 

s s « 

O’ Q 
ssi m >t 
H O Q< 
CD n o 
P« » H 
OOP 

0 o 

1 ^ I 

sis. 

3 § S 

>o o 

i a ■§ 

OHM 
Q O O 
Oi XD a 

§!•§ 

P o *4 


P P 

i a a 
a as 

o o 

S H « 
OOP 

0 

« « tt 
o d 
s> P p 

a S- 

a o“ 

p ■ 

S ^ 

P ts 3 O 
CD CO p 

H « >* 

1 1^ 

« 3* o 


p o g 6 
, P 0 o 
HOP 
P o cq 
_ « P< o 
o a H 
o a- H p 

P -o O H 
P O ?> O 


9 o 

POP 

0 n 

1 i ^ 

5 p o 

OOP 
H n o 

s a a 

p ^ ^ 

3 ii 

H 2 p 

I o p p 

O tc 

Q p 


p a 

£! a 

a o ^ 

2 s?? 

p ^o 
O P 


^ O tA 
P O 
P« o 
a o 
p« p a 


a s a 

O P 

pop 

CD H ^ 

s a a 

0 & § 

H 0) P 

1 s ^ 

OOP 

p a p 

S a 
t>* o 


5 i 


■p 

1 ^ 

o 

o 

dc 

Q 

e 

•H 

5 ? 

f 

j 2 

fco 

o 

o 

*d 


Q 

C3 

-P *d 

S 

A 

rH 

O 

•o 

4-> 

a 

4^ 

<Q 

O 

*d 

H 

tc S 

3 

i 

2 

o 

rS 

a 

d 

6 

o 

9 

> 

-a 

M 

o 

A 

& 

\ 

o 

n 

o 

P 

*ri 

o 

-p 

ca 

3 

-H 

to 

■a i 

w 

Q 

i 


rH 

o, 

d 

§• 

o 

CD 



1 

9 

•H 

O 

9 

9 

to 


3 o 

d 


P5(3B jg 

Of' POOR QUALITY 



A 

<D 

S SvJ•;k^ 

i:i M w § 
o t-* « 

03 33 U ^ 'N 

* g Sl^'^vTN 



w| 

ads 


\<j , 

pCi ^<1> *W^ 

o o S >c C M 

g g 

1 1 a 

la 'V' 5 ?^t 

O ^ viw t ^ PJ 

§■ I 1 

_ a I 


Vi cA *m 

w >o o 

Q) W I j ■ I 

1 1 

o w *3 t S> K 

« « Q rii M 


Ij::^ ‘ I 

vVI *k _ 



I 


' i 5 ! ^ 

^ ** , I p- i— ^ 

g II 


t:u 


li4 

SiT? 



- ^s; s 

^ gl l rcs4*'^ii 


a 

» 


Sci?g«f 


o 

» 

CD 


M W 
It 





6 ^ 

H O 

o o« 

M CM rt CO*CM cnti ^ 

p« . * . . ^ 

¥ *ha -f* ¥ +P*^ » 

X a 
« * 

4:^*4 -4'. 

, ,CM CM * (t) ^ 



4 o 


OH H 


H M H 


H H 


C^ 





» &NJ 


03 & a 

I § 

03 3 9 “ 

65 *d 

jq o ^ o CQ 

^ g£ c 3 '-' 

« g"S 

a w K 

S " S 5 

n . « 3 


pa< 


CM 


•f 


^ k * fc. ***55 

S g^cS t 

O w C^O ^ 
(O Ml) b 5 

S B « § 


a 


a s< 

O Cl 
CO o 


d o H 

" E a 

o 


J i«p 


1.^3 

IS S 



g 


d 


r» d S 

lo’ S S 03 S 

OH H 
CO ID 

§ & a 
<o H d >o 

_ p( ttf d p 

"a ip 

O H P* 

P 4 U P 


i 098 


I 



H 

1 

1 § s 

frt fri 

CD 

B 8 

i 

o 

CJ **♦ 

n •< »4 
5? r*^ n 

^2i 

B^l 

O O 
H >0 

S.O 

1 







i* 

H 3 

& s 




I 

I 


u 

o 

« 

o 

g 

s 

o 

a, 

M 

N 

H 

o 

o 

n 

a 

a 


M 

o 

o 

CD 

Vt 


H CD 

o e 

(D » 

Si 

o O 


X H ><^ 

fXJ 


a « . 
s ' 

■Lo5: . 

o -5>l' 

’>«£. ' 5 . 


X 


Jl ^11 


'■^ V5> _ » ^ ® 


a ^ 

9 s 

n o 

° 1 
o 

S m 

8 B 



I |t fl C?J p 


1 ^ » 
a » . 1 . 


u> n 
\ It 

‘3 is *'■ r S 

C 9 t 4 I O 

:ri 

& a a ^ gt 

o g g w *s^«s 

t? 34 <D CO ^ o 

S a ^ K \ a 

o N cp o A 

. ^ M O H H : H 

^ ^ n d O n <0 

kj I-- § g ?. 5 

l&lls.t <l)g 

=■ § “ «?«1 

o o » 0*0 ^U4 

a«apisi (4 *<y^ 

s I s s i & I 3 

SAsag -uNS 

al?e< a ^ 3 


a^g g=i i I 

a a a • s 

g'o O 9 g 0 - 

L rgv5= a ^ 

8 A.S « s 

C« O P* 7 ^ d d 

34 §|i"; “|i 

o a o fcf o 


kp\]S 




w ^ Vc 

*3 ^ >r 


f 4 'O 


° < 

S>NJ 


<! . 
^S 5 


<J§ 

- "I 

.‘?i!-i «:'! 


iW| 




^3 ^8 


^4 


tiS II ^ 


H C 9 

« g S, 

^|s 

^ 

o d <d 
o d 

a d o 
odd 
M 14 d 
dud 

9 ^ & 
3 g 

I tt; 3 H 


3 d o 

^ d 3 

o fri Pi 


no* 

3 P- 

M *»!' 

U H P 
C 3 O p • 

3 a aa 

w Si a ’‘^ 


“f- <' 

*"iWS 

'ij 




CS'ilMJiWJCb’ , ^ 

3 ." II ^11 -^" li^ 




a 

OJ*^ ^ 

ra ^ ^ 


CO O §4 
K O O 




d o o 
H :s 5 

Odd 
3 0 0 
H O* ?=J 


o d o 

^ P f 4 
Odd 

5 a 

0 cd o, 
d d a 
H O 

1 I a 

d ^ s 

d o o 
« d 


d d 

<□ o d 

cl CJ ^ 


:yy/ 

3 ^ 

« Ll^ ^ 


o d d 

w a Q o 

5 1 *}^ 1 1 

cd a w * 01/ 

d *s _ CO 

I 

g „ a 


i^g^, 

d ^ a 
o Ic 

S*8i ^ 

, Eh Cnrv<_, 

Oi 

iS S'-- 


CiV a 
1 


R ^ n. 

o 

>< as 

*c O 'S 

> o n 
V 3 I 

V Sfuj 

S 5 C 


v-^N 

i-^v 


i I 


V « 1 

csi' ^ w 

^ *s. a M 

B M 

I|’ <0 ^ 

I* ;a 


ai 

o 





^ o O 

§ “ I 


g ^ , 

goo 

i i 1 3 


" o 
•o fcJ 

g . 

o 31 
V> 

^.58 a 
s. g 
s a 

2 p>> 3 


i*< 

II I 

? " i V 


| 3 g:- 

O H H ff 5 * 

K a “'T' JO 


I 

ll;l: 

_ o Lj n 
^ O w 

i IP} “ 3 


*A S 3 

H o XV 
o a i 

I 


a o 

I ^ 

CJ o 

e a 

in i 3 

^ a & 

Q « 

J! 5 ^ & 


''‘^'N O 3 sj « lO >> Si 

v5^ I a> 

|sqp-ivi»| 


o • s-^ 


» I CS' s. 

i ^ a Jim § 

. 0 i^ s g 

»i 3 -i S JJ 

« •Op 

O ti B5 •» 

& S S H 

g 8 Ci 

^ K K H 

i{ H <3 3 


S 


5 '■'g 


I i« - I 
s.'&i ia.il 




CJ w g 

o g ^ 

""S3 


- y 


W ii 5 < 

i ItM 

'o* ll-" 

- S ^ 


- Td & s a 

oro |5 «•*-! ^ - 

3 S i<r~t 

” 3 m • <o *>I^ 


^ a 

2 g 

s a 


<5? 

^Nf 

V 

n 

i ^ 

s<, 


o 31 « 

S SO 

CO a Q 


^ « 

^ S K 


V: 


3** »«• 

a<vO 
^ ( 
Sn 

Q ^ 


'"n i.^r 
^L|' 3^ 

• Cfl 
^ t-i 


i ^ I fiC re? i 

g,-cA Ci/ g S 

« ^ - 9 a 


0 V 

H ^ S 

1 <:s»« S 

S fs' 2 


S « » 


CO H g 
Pi Id » 

^ o o 


a 


n ^ <i 


S S’ a >5 
K P< H g 
o o X 
- 3t 

<D 1 H * 

a O 

rr 

g cl • 

3 ** 9 

a rc 

grO « s 
o * xj 3 

® PUli ? 
R F< ^ n 


,^«i 


V 

<r 

<i£^ 

.>tsj> 

^ 4 - ** 
cqj 

^K| 

I '• 

« Jr' 

vel;^ ^ 

II ' 


3 to 

.— » H 

s 

JU3 

- 


O If a *‘ 

lv 

"J I*'.C 

c 3 ^ o 

p='==y V 

o my^ 

a aw* 


1^53 V 

at 


,.k m 

^ - 
Oi 

’^wr. 


r-i ^ 


g'Xi 


.Sf V| 
'^'<1 § 
<a ^ « 
H ^ s 
< g “ 

V . 9 M 


^«S S + 
+ 

^ ti a H 

. p 

2 H 

‘"WSS“ 

II a s 


>:DsJ> 

*v- 


a 54 


I 

*Vd 

'S)^ It 

O *o 

<1 ti 


i 

‘ 


Iv2- 

a uo 

a ^ 




O^ 

« ^ 
o 




« 2 

•s ** 

N g 

3 I 

3 ^ 

la 


t 

1 


H M 
O 

w 

to 


u: 


§ 3 

&i 

O rt 

3 s 

g ^ 

§ » 

. Q o 

tt CJ 
-< !-• 

« O 
tr H 
O tr 


M 

O 


Vvj 

^>0 

ll 

r-» 

<=} 

II 

r-i 

S 


hT 

It 

r~T 

U 



5 3 »§*^ 

OH 


^ ^ 3 S S 

CO ^ t ^ H S 

o* CD cn 0 


^ a ^ p 


2 ? 
o SI 



q g 

■A I ^ B 

M pi 




I A ttj , ^ 

^ « g A g s 

§ . a & g s 

8S 3 as 1^ 
p, a w 

o o R iJ 

a 5 o a 

«agi ■ 
iliSrs, 


» 

II. 

s a 

O 


S 

H 



It® il.i 9 

w ^ g ft 

o & 


o 

s 

a 

«} 

o 

g 


w * H* 
® Cj 


in , 

as 


« s _ p,, 

§ a ^ >a g “ a ' o 

gas-”,. ° a 

Aoob*’* 

wosa irS*- 

« a H S « 

IS| CJ * P H O 

O g o g 

0 14 c;l 

OKU 

1 5 S 

OHO 


,r-nO 


ft B S « s 

B a H 

^ § E 8 

o g B ^ 

« a a m _ 

H H O 


w 

s s s 


i' 1 '^ 



raa 

g^-i s 

a S' 

si « 

a g 

& s ^ 

O 0 

o N o 
3 . ^ 

o' o a 

« w 5 * 





a& 

63 R 

^ s 


‘M 


(« 

o 

S 3 ’ 2 

s & 

o o 


cQ m 


t 


